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COLLISION INTEGRALS FOR NONELASTIC PROCESSES 


IN PLASMA KINETIC THEORY 
By John R. Viegas 
Ames Research Center 


I . SUMMARY 


Often in the treatment of plasma dynamic problems the effects of 
nonelastic encounters must be included. In this note collision integrals for 
nonelastic processes that are important in plasma kinetic theory are developed. 
The following nonelastic interactions associated with the electron Boltzmann 
equation are presented in detail: 

(1) Inelastic and superelastic encounters between electrons and neutral 
species , 

(2) Ionization and three-body recombination encounters in which the 
third body is an electron, and 

(3) Photoionization and two-body recombination encounters. 

In addition, collision integrals are given for some of the processes that 
affect the transfer equations of heavy particles (ions and neutrals) and 
photons. Included here are radiative excitation and de-excitation encounters. 
The compatibility of these collision integrals with the forms associated with 
the usual macroscopic rate equations are demonstrated. 


II. INTRODUCTION 


Often in the treatment of plasma dynamic problems the effects of 
nonelastic encounters must be included. Such effects are frequently of primary 
importance in the calculation of electron distribution functions and in the 
evaluation of electron continuity and energy equations (refs. 1-13). To a 
lesser extent they can also play a role in transport property calculations 
(refs. 1, 7, 13). 

The purpose of this note is to present, under one report, in a consistent 
form, collision integrals for nonelastic processes that are important in 
plasma gasdynamic situations. These collision integrals will be presented in 
the format usually associated with elastic collision integrals (ref. 14) . It 
is hoped that this format will lend itself to quick assimilation by the reader. 

This paper represents some modifications and extensions of the work of 
many authors. In addition to some of the previously cited references, the 
work of Fowler (ref. 15) was found to be particularly helpful. Some of the 



material presented in Section IV has also been published in reference 16. 
However, for convenience, it seemed appropriate to reproduce the material here. 

This paper is organized as follows: The first part presents generalized 

nonelastic collision integrals for electrons that account for 

(1) Inelastic encounters between electrons and neutral species, 

(2) Ionization and three-body recombination encounters in which the 
third body is an electron, and 

(3) Photoionization and two-body recombination encounters. 

Next the nonelastic collision integrals that affect the distribution of 
neutral species, ionic species, and photons are treated briefly. Included are 
radiative excitation and de-excitation encounters. In the last section the 
compatibility of these collision integrals with the forms associated with the 
macroscopic rate equations is demonstrated. 


III. COLLISION INTEGRALS FOR NONELASTIC PROCESSES 
ASSOCIATED WITH THE ELECTRON BOLTZMANN EQUATION 


In a collision-dominated plasma composed of electrons, neutrals, and ions, 
the rate of change of the electron distribution function^ fCvs) as a result 
of nonelastic encounters can be represented as 


/^e^s\ _ /^e^s\ /^e^s\ /^e^s\ 

V /ne ^ V 'Wex V /ion " V ; 


Ph 


( 1 ) 


where f^ = f(Vg). The three collision terms on the right-hand side of 
equation (1) represent the nonelastic interactions considered most important 
in collision-dominated plasmas . They, respectively, account for 


(1) Inelastic (excitation) and superelastic (de-excitation) encounters 
between electrons and neutral species. 


(2) Ionization and three-body recombination encounters in which the 
third body is an electron, and 


(3) Photoionization and two-body recombination encounters. 


In this section the collision integrals for these terms will be developed. 

The state of the heavy particles will be indicated by the subscript j, k, or 


^As is common practice in the treatment of collisions, the spatial and 
temporal dependence of the distribution functions are not shown explicitly. 
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Z added to the species subscript. These states are considered to be 
specified by the principal and total angular-momentum quantum numbers and 
degeneracies. 


Inelastic and Superelastic Collisions 

In this section collision integrals are developed for the conservation of 
electrons in the set v^, ^ result of inelastic and superelastic 

encounters that cause neutrals to undergo excitation and de-excitation reac- 
tions, respectively. A typical interaction can be represented by 

1 Ic 

e + (Atom)^ e + (Atom)^ (2) 

where the subscript n stands for the type of neutral atoms and the super- 
scripts j and k (k > j) correspond to the state of the neutral particle. 

This type of encounter can be distinguished by whether the reaction 
causes the electron to be subtracted from or added to the set. When an excita- 
tion (j -> k) of a neutral atom is initiated by an electron with velocity v^, 
the energy equation is 

(l/2)nigV52 + (l/2)mnVnj^ = (l/2)nigV^2 + (i/2)mnVnj^^ + (3) 

Here Vq is the electron velocity after the collision, Vj^. and Vn^ 

j ^ 

velocities of the neutral particle before and after the encounter, respectively 
(the subscripts j and k serve the dual role of distinguishing between the 
neutral before and after collision and representing its state), and A., e hv., 

r, . 

is the excitation energy for the interaction. When an excitation results in 
an electron being added to the set, the energy equation is 

(l/2)mgVi2 + (l/2)mjjVj^^^ = (l/2)mgV^2 + (i/2)mnVj^^2 + (4) 

where is the velocity of the electron before the inelastic collision. 

The inverse (superelastic) encounters, which also contribute to the number of 
electrons in the set, are described energetically by equations (3) and (4). 

The number of electrons lost to the set Vg, d^Vg in unit time per unit 
volume as a result of inelastic collisions with neutrals of velocity range 
Vnj> d^Vrij , such that the neutrals are placed in the set Vnj^, and the 

electrons are scattered through an angle Xqs> measured relative to 

Vg - Vnj in the center-of-mass frame, with the velocity Vq, d^Vg is 


3 



fs ^\^n. ^Vgsj^j 


( 5 ) 


This type of transition is written: 


\ 


In the preceding inelastic loss expression, ^5 j CSg j 


is the differ- 


ential scattering cross section for the excitation collision in question. 

Iv^ “ J , fi - 


The 


relative velocity g^j is defined by g^^ ~ 

velocity distribution function for the species n 
normalized on the number density for this state. 


n, = ^nCvnJ is the 


‘3 ■ “J “ 
neutral in the state 


In writing g^ 


a . , we have addressed ourselves 


^ and Xos arguments of 
s j 

to an examination of excitation collisions which cause heavy particle state 
changes that can be classified by the principal and total angular-momentum 
quantum numbers and degeneracies . The cross section should have enough infor- 
mation specified in its arguments to determine, along with the conservation 
equations, the velocities of both particles in the appropriate reference frame 
after the collision. In the dynamics of inelastic collisions in the center- 
of-mass coordinates, changes in the component contributions to the angular- 
momentum vector of the system will cause the plane of the particles after the 
encounter to differ from the plane of the particles before the encounter. If 
no strong external magnetic fields exist, the direction of the total orbital 
angular-momentum vector of the heavy particles can be taken to be arbitrary 
both before and after the collision. Only the magnitude of this vector is of 
interest. Then by symmetry the azimuthal dependence of the orientation of the 
relative velocity vector after the collision with respect to g^j can 

be neglected. Thus only the angle Xqs between and g^j needs to be 

specified to describe adequately our inelastic collisions. 

This argument can also be illustrated by the interaction between two 
monoenergetic streams, one of which is taken to be the set of electrons having 
the velocity range v^, other is a stream of neutrals in the set 

Vnj, a result of the averaging effects of the streams the sum of the 

total angular- momentum vectors of the heavy particles can be expected to be 
zero. This initial symmetry would be preserved after the encounter allowing 
us to neglect the azimuthal dependence of g^j^ relative to gg j • Thus, only 

the single deflection angle Xqs needed as a parameter of the collision. 
Knowing Xos^g • using symmetry considerations one can find the velocity 

of the particles after the collision in the center-of-mass reference frame. 

These same arguments also apply for the other inelastic (superelastic) 
collisions discussed in this paper. 
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The number o£ electrons lost to the set Vg, (Tv^ in unit time per unit 
volume as a result of superelastic collisions described by ^ Vi 

is 


1 


U.1LJ. U 




( 6 ) 


Here (g , ,x is the differential cross section for this super- 

SK SK is is 

elastic encounter. The angles are measured in the center-of-mass reference 
frame relative to number of electrons lost to the set per unit 

time per unit volume is then the sum of (5) and (6) . 


Similarly, the number of electrons gained by the set per unit volume per 
unit time as a result of the inverse encounters represented by 
W’X (v^,Vn Vg.Vn J is 


and 


f O d3 (g^j^ , X^3 ) dfi^^ 


f 1 d3v,fn. d^v^.g^jO^ . Cg^ . ,Xi3)df^^^ 


(7) 


C8)- 


Expression (7) is the (superelastic) inverse of (5) . Expression (8) represents 
the inverse of (6). Typically the relative velocity is defined by 

§ok " - %!• 

A simplified expression is obtained for the net gain of electrons to the 
set per unit volume for the reaction given by (2) if the principle of detailed 
balancing is applied to the inverses among expressions (5), (6), (7), and (8) 
before the integrations are made over scattering angles and particle velocities. 
Applying this principle, we equate (5) and its inverse (7) at thermodynamic 
equilibrium to obtain 


The equilibrium distribution function for fg is given as 

/ me \3/2 -mev2/2kT 

(^s)Eq (^e)Eq l2TrkT) ® 


(9) 


( 10 ) 


^To be consistent these angles should be labeled Xe , dfi and x , dS^ • 

so SI SI 

however, it is easy to deduce that for inverse collisions as described here 

X = X 9-nd X = X . 

'^so ^os ^si ^is 
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where T is the equilibrium temperature, with corresponding expressions for 
fo. fnk« and • • The relative populations o£ various excited levels of an 

atom at equilibrium are given by the Boltzmann distribution 



(11) 


where and ojj are the atomic degeneracies associated with the states 

j and k of the neutral particles. The differential velocity elements in 
equation (9) can be related via their Jacobian as 


d^Vnj^d^Vo = ^|^^d3vnj 

Then the combination of equations (10), (11), (12), and (3) with (9) 
following detailed balancing result: 


( 12 ) 


yields the 


03 . ^ 

^ok^ok ^^ok^’^os^ ” o)j^ ^sj^sj ^^sj '’^os^ 


(13) 


A similar detailed balancing analysis with (6) and (8) yields 


^sk^sk ^^sk^^is^ ” 0), j^l j j 


(14) 


The differential velocity elements for this encounter are related by 


(15) 

If the electron velocities before a superelastic and an inelastic 
encounter are denoted by Vg and Vj, respectively, it is easy to show that 

equations (13) and (14), relating the differential cross sections, can both be 
written as the following general detailed balancing result: 


g?-i°T-i (St.-’X) 




The terms ggj^ and g^^ in equation (16) are defined as 


>Sk 


= - Vn,, 


and 


( 16 ) 


6 



glj = IVj - Vn. 


The velocity elements in this general notation are related by 




^Sk 


d^Vj^. d^Vj 


( 17 ) 


Equations (12), (13), (14), and (15) can be combined with (5), (6), (7), 
and (8) to yield an expression for the net number of electrons gained by the 
set Vg, d^Vg by the reactions represented by . ) and 

(^1 ,Vn. Vs,Vy^j^) . When the resulting expression is then integrated over all 

possible scattering angles and heavy particle velocities, the following expres- 
sion is obtained for the net gain of electrons to the set per unit volume per 
unit time for the reaction given by expression (2) : 



j^k 


d3vc 


r//(' 


— — fs^n )g .a . d^v 
i o Hr ( 0 ^ ^ "jpsj sj os rij 




CIS) 


The sum of (18) over the various neutral species and states of these neutral 
species divided by d^Vs is 



(19) 


In this equation the first collision operator represents the superelastic gain 
and the inelastic loss to the electron set v^, d^v^ for all possible species 

n and for all possible values of j and k. The second collision operator 
represents the inelastic gain and the superelastic loss to this set for all 
possible species n and for all possible values of j and k. 
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Ionization and Three-Body Recombination Encounters 


In this section collision integrals are developed for the conservation of 
electrons in the set during collisional ionization and three-body 

recombination encounters. A typical reaction in this case can be represented 
by 

e + (Atom)^ e + e + (lon)^ (20) 


Here, as previously, the subscript n corresponds to the type of neutral atom, 
and the superscripts j and I signify the state of the neutral atom and its 
related ion, respectively. 


When the way that electrons enter or leave the set Vg, d^Vg is of 
concern, the energy equations for this type of encounter become: 


r2 = 1 


2 "'eVs" + 2 = 2 ’"eV2 + 2 ^ 2 

when the ionization is caused by a Vg electron and 

1 2I2I 2I 2I2. 

2 "'eVi 2 "^nVn. = 2 + 2 "'eVs + 2 ^ 


( 21 ) 


( 22 ) 


when the ionization results in an electron in the set. In these energy equa- 
tions nij^ and are the mass and velocity, respectively, of the ionized 

type n neutral atom in the state distinguished by the subscript £, and 
A. is the ionization potential of the atom from its jth state of excitation 

to its related ion in the £th state. In equation (21) V 2 and v^ are the 
velocities of the two electrons resulting from the ionization. Making use of 
the principle of indistinguishability (ref. 15) we will not distinguish 
between bound and unbound electrons during an ionization or three-body recom- 
bination encounter. Thus, the electron that corresponded to the ionizing 
electron after the ionization will not be specified; it will merely be stated 
that two electrons result from the ionization. Likewise in a recombination 
encounter, no attempt is made to identify which of the two free electrons 
becomes bound and which remains free. Similarly, in equation (22) a Vj 
electron caused an ionization which resulted in two electrons, one with veloc- 
ity V 3 and the other with velocity V 3 . It should be noted that the V| 
appearing in this section is not the vx that appears in Section III (p. 3). 
They are distinguished by the type interactions being considered. 

The number of electrons lost to the set V 3 d^V 3 in unit time per unit 
volume as a result of ionization encounters with neutrals, such that the 
following transition occurs 



^ V2,V^,Vi^) 
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is 






vn.g3, 


£ / 
.a . (| 




42 


CM . 


42cM 2S 


( 23 ) 


Here 


"sj'SsjSXzs 




is defined as the differential "cross 


section" for an ionization encounter between electrons in the set 
and neutrals in the set v, 


Vs, d^Vg 


ny resulting in two electrons, their veloc- 

ities given by V2, d^V2 and V4, d^V4, and an ion in the state £ with its 
velocity range and state signified by d^Vi^. Here the semicolon in the 




argument of the cross section separates before collision parameters from after 
collision parameters. The angle Xoc is measured relative to the direction 
of g^j in the center-of-mass frame. The velocity *'^ 42 ^j^ is defined by 

^ where X42^j^ is direction of V4 measured relative to V2 

in the center-of-mass frame. The center-of-mass subscript CM is used here 
to avoid confusion with non-CM velocities. 


Again symmetry arguments have been used to arrive at the above choice of 
cross-section parameters. As presented, sufficient information has been 
specified that the velocities of the particles after the encounter can be 
determined when symmetry considerations are included with the conservation 
equations and the initial velocities. There is no single "correct" way to 
designate a differential cross section of this complexity. For example, let 
us examine qualitatively the dynamics of an ionizing collision as regards a 
choice of parameters for the cross section. Consider the interaction of a 
monoenergetic beam of electrons with a monoenergetic beam of neutrals which 
results in ionization of the neutrals. Before the ionization encounter 
between an electron and a neutral particle, the two particles lie in a single 
plane in their center-of-mass reference frame. After the encounter, however, 
the resulting three particles are not restricted to move in a single plane in 
this reference frame. Averaging all such encounters in the beams will still 
give a symmetry of sorts about the relative velocity vector that is, the 

resultant momentum vector of any two of the scattered particles (say the elec- 
trons) from a single collision can be found equally likely in any azimuth about 
ggj as will the momentum vector of the third particle. These momentum vec- 
tors, that of the above resultant and that of the third particle, are equal, 
antiparallel, and lie in a plane which has the same symmetry about g • as 

the particles on page 4 after an excitation collision. This symmetry can be 
represented here if any of the three velocity vectors is allowed to be arbi- 




We choose v- 


trarily located in an axial sense with respect to g 
to have this symmetric character and thus only specify X2S 


£ 

y . , 

SI 


-CM 


The 


collision is then completely described by specifying the velocity of one of 
the remaining particles relative to We picked V4 for this distinc- 


tion. It is important to note the arbitrariness of these choices. Thus some 
other choice of parameters could have been made. It may appear from the above 
choice that the "created" electron must have the velocity V4. This is not 
the case. Consistent with the principle of indistinguishability we have 


9 


merely stated the velocity of one of the resulting particles without 
specifying which particle it was prior to the encounter. 


The number of electrons lost to the set 


->■ 

Vs. 




Q Vg in unit time per unit 


volume as a result of three-body recombination encounters with ions and other 
electrons such that the v^^v^^^) transition occurs is (ref. 15) 








gog 


3*S 




I dS7 


]s 




The relative velocities are defined as 






and 




1^3 - Vi^l 


Two cross sections appear in expression (24) because there are two ways in 
which electrons can be lost to the set d^Vg by the (Vg.v^^vj^ -> v^^v^O 

reaction. That is, on momentarily distinguishing between electrons, either 

the Vg electron can become bound to the ion and the V 3 electron becomes 

the vi electron or the Vg electron can become the v^ electron while the 

Vo electron becomes bound to the ion. Thus o3 is the ^^differential 

^ S35, is . 

cross section^' for the recombination encounter in which s 1 and dQ, ^ 

3s£ 1 3 

is the ^’differential cross section” for the recombination encounter in which 
3 1. From the principle of indistinguishability we can say 


S3£ 


dn 


IS 


= 


3S£ 


dQ 
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Thus, expression (24) becomes 


We 


2f3d3v3f3d3v3fijd3vi^g^^gj^<,3^^(g^^,g3j,Jg^g^;Xis)dS,s(l * C25) 

have taken „(g ,g^o>X as the "differential cross section" 

S 3 £\*S£’* 3 £’^g^g^’^lS/ Is 

for this recombination collision. After this encounter the remaining free 

electron is scattered into the angle x, measured relative to 

^ 6 ^is’ is sc^^ 

The angle between g^^ and g^^ is given by Xg g • Again, the choice of 

3 S 

argument for the recombination cross section is not unique, and also we do not 
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specify which of the two free electrons becomes bound (the cross section is 
symmetric in the subscripts s and 3) . The delta function 6-v •> is included 

in (25) to account for the case when both electrons participating in the 
recombination encounter are in the set ^g, d^Vg. This function is defined as 
follows : 


|o ; 

’ ”3 

Thus, for any smooth finite function K(^ 3 ) 


6^ ->- = 
V3,Vs 




V, 


Vc 






d^Vg 


0 


so that for practical purposes the delta function can be omitted. In an 
analogous manner, the number of electrons gained by the set per unit volume 
per unit time as a result of the inverse recombination (v 2 ,v^^,Vi^ ^ Vs^Vj^^.) 
and ionization ^ Vg,V 3 ,v^^) interactions are, respectively, ^ 


fg dSv^f^ d^V^fi^ ;Xs2)d^^S2 


and 


2fi d3v,f„. 


dQ^, d^v,,^ (l + 6^ 


CM - ^35 


fl + 6 ^ (27) 

CmV 


In expression (27) the coefficient 2 accounts for the two ways in which elec- 
trons can enter the set via the (v^,Vj^. interaction, that is, 

either ^ v^ and the bound electron becomes ^ 3 , or ^ V 3 and the bound 
electron becomes V 3 , and the cross sections for these reactions, which are 
symmetric in the subscripts s and 3, have been set equal on the basis of 
indistinguishability arguments. As a result of the collisions being inverses, 
and the symmetry discussion following expression (23), the scattering angles 
in (26) and (27) are the same as those for the direct collisions and will sub- 
sequently be treated accordingly. Note that the angles are measured relative 
to velocity sets and not relative to particular particles. The delta function 
in (27) is to account for the possibility of the ionization resulting in two 
electrons entering the set. 

As in the previous section we will derive a simplified expression for the 
net gain of electrons to the set per unit volume for the reaction given by 
(20) through the principle of detailed balancing. Applying this principle 
expression (23) and its inverse (26) are equated at equilibrium to obtain 
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At equilibrium the distribution functions are Maxwellian: 




= (ne) 


Eq 




m = s, 2,4 


(29) 


where g = m3/2kT, with corresponding expressions for fj^. and f^ . The 

3 ^ 

number densities of the ions, the neutrals and the electrons are related at 
equilibrium via the Saha equation, which can be written 



me ' 


3/2 


( 30)3 


In equation (30) coj and o)^ are the degeneracies for the type n atom in its 

jth excitation level and its related ion in the £th excitation level. Com- 
bining equations (21), (29), and (30) with (28) the following detailed 
balancing result is obtained: 


d3vn- 


(31) 


where 




A similar detailed balancing analysis with expressions (25) and (27) 
yields 


= Hg^ .ajj d3v33^^ d3y^ dSy^ (32) 

As for inelastic and superelastic encounters, equations (31) and (32) can 
both be included in a single expression. If subscript Z denotes the elec- 
tron causing the ionization and subscripts R and B denote the electrons 
resulting from the ionization, equations (31) and (32) can, consistent with 
previous notation, be included in the following equation: 


Since the following relations between the differential velocity elements hold: 

d\d3vg(#Vi^ = d3cd3gj^^d3gg^ 

3as written equation (30) includes the approximation that m^ = mj^. 
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and 


d^Vnj = d^Gd^g^j 

where S is the velocity of the center of mass, equation (33) can be written 
as 






CM 


Zj 


(34) 


Equations (31) and (32) can be combined with (23), (25), (26), and (27) 
to yield an expression for the net number of electrons gained by the set 

Vg, d^Vg as a result of the interactions represented by (Vg,Vn- "^^2’^4’^i ^ 
* D ^ 

and (v ,v ,Vn.). This resulting expression can then be integrated 

over all possible scattering angles and heavy particle velocities to obtain 





(35) 


This equation represents the net gain of electrons to the set per unit volume 
per unit time for the reaction given by expression (20) . 

Now, on summing over the relevant neutral species and the excited states 
of the neutral species and its ion and dividing by d^v^, we obtain the follow- 
ing collision integral for the net gain of electrons to the set per unit 
volume of phase space per unit time as a result of collisional ionization and 
three-body recombination encounters: 


Peh\ 

V Aon 


E 

n,Jl, j 


[ffA 




n£ 


H 


f^f. 


nOssj 


SJ 


d^^ d^v 


2S 


42 


CM 


d^Vr 


Here the first collision integral represents the gain by three-body recombina- 
tion and the loss by ionization to the set v^, d^v^ ; the second collision 
integral represents the gain by ionization and the loss by three-body 
recombination to the set. 
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Photoionization and Two-Body Recombination Encounters 

The remaining nonelastic collision term is that associated with 
photoionization. Physically, a photon of energy hv is absorbed by an atom 
in the state j, resulting in an ion in the state £ and a free electron. 

The reverse reaction is when a free electron and an ion combine with the 
emission of radiation. A typical reaction for this case can be represented by 

e + (lon)^ (Atom)^ + photon (37) 

In terms of the relative velocity g^^ the energy equation for this 
encounter is 


hv 


1 2 
2 ^i^s£ 


+ A. 


(38) 


where is the reduced mass defined by m^m^^/ (m^ + m^) . We consider only 

nonrelativistic electrons in this analysis and neglect the momentum of the 
photons relative to the momentum of the electrons or the heavy particles. 

Then in the center-of-mass reference frame the direct encounter (photoioniza- 
tion) will result in an electron and an ion moving in an antiparallel direc- 
tion arbitrarily oriented relative to the direction of the incoming photon. 

If we consider the interaction between a monoenergetic photon beam and a cloud 
of neutrals in the state j in the center-of-mass frame which are not polar- 
ized by an external magnetic field, then by symmetry arguments we can conclude 
that the collision probabilities will be independent of the ''scattering’^ angle 
of the recoiling particles. Thus, in any direct collision v would be the 
only parameter needed, in addition to the conservation equations, to determine 
the velocities of the resulting ion and electron in a particle direction. 


Let Qj(v)df ^5 represent the differential cross section for the ioniza- 
tion of a neutral atom in state j and set Vj^. by the absorption of 


radiation of frequency v resulting in an electron being emitted into the 
angle measured from some arbitrary reference axis with the speed Vg, dVg. 

The resulting ion will be in the state £ with velocity Vj^ , d^Vj^ . The 

solid angle associated with the ion motion will be in the center-of-mass 

frame. The number of electrons gained by the set V 5 , d^v^ per unit time per 
unit volume as a result of photoionization can then be represented by 


fnj d3v^^f^(v)dvcQj dfig 


(39) 


In (39) is the photon distribution function per unit frequency per unit 
volume; it is normalized on the photon number density and is related to T(v), 
the specific intensity of v radiation integrated over all solid angles, by 
(ref. 17) 




I(V) 

chv 


(40) 
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T(v) is proportional to the radiant energy density at this frequency. Here, 
and in the previous expression, c is the speed of light. It plays the role 
of the relative speed in expression (39). Note that we could have included 
explicitly here the direction of propagation of photons Xv worked with 

f^C'^ 5 X^)dv rather than 7^(v)dv. This will be done in the discussion on 

radiative excitation and de-excitation in the following section. These 
distribution functions are related by 

•'o 


where d^^^ is the solid angle centered about x^- 


Two-body recombination, the inverse of the above encounter, results in a 
depletion of electrons from the set. An examination of the dynamical equa- 
tions in the center-of-mass frame will reveal that the only parameter needed 
to determine the frequency of the radiant energy is the relative speed 

Then, defining P^Cgcp) the differential cross section for radiative cap- 

X, ^ ^ 


ture (two-body recombination) of an electron in the velocity range Vc 
by an ion in the set 


d^v. 


'S^ ’'S 

d^Vj^^ which results in the emission of a quantum 
of radiation hv and a neutral atom in the set ^nj ^ expression for 

the number of electrons lost by radiative capture for the reaction represented 
by (37) can be written as 


d^ 


Vcfi 
S Ic 


d^A 


,gsJiPo Cg 


sZ ■ 


(41) 


Electron captures result in both spontaneous and stimulated emission. 
Thus can be written 


pj(gsj^) = 4^(g3^) H- I(v)3j(gsj^) 


(42) 


where the first term on the right-hand side is the contribution from sponta- 
neous capture and the remaining term is the contribution from induced capture. 


In order to obtain a simplified expression for the collision integral 
associated with (37) we apply the principle of detailed balancing and equate 
(39) and (41) at equilibrium. Utilizing (42) and (40) with this equality we 
obtain 






4ira^ + 


[I(v)] 


Eq 


d3 


Vc 




[fnjI(v)]Eq 

hv 


Qf(v)d3vn. dv 


(43) 


At equilibrium the distribution functions are Maxwellian and are given by (29). 
The equilibrium specific radiation intensity, given below, is Planck’s black 
body intensity: 
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TTr n _ STrhv^ / 1 \ 

= ^2 (^hvAT _ J 


( 44 ) 


The Saha equation (30) relates the number densities in equation (43). Com- 
bining equations (29) , (30) , and (44) with (43) , we obtain 

e'^^/’''^ejjd3v3 d3vi^ = qJ ( v) ^ d3vn. dv dfi 


Since this equation must be independent of the temperature, the following 
relations must hold: 

i - . _ 2hv3 j 

and 

qj(v)dfi 3 dv d3vj,. = H'^ ^ g^^aJ(g3j^)d3vg d3vi^ 

The differential elements in equation (47) can be related through the 
Jacobian of the following transformation: 

dn^ dv d3vn. = |j|d3vc d3vj 


s 

(45) 


(46) 


(47) 


1 


(48) 


To evaluate |j| we make use of the following differential relationship, which 
holds for our photoionizing model: 


When this relation is combined with (48) the result can be reduced to 

dag dv = |j|d3g3^ (49) 

Since d3gg^ = g|^ ^^s’ written as 


dv = 


dg 


sZ 


so that 


J = 


9v 


®S£ ^®S£ 


Then using (38) we can find the following relationship 
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8 V 


h 


Thus the Jacobian becomes 


and (48) can then be written 


J| 




dfi dv d^Vr, 

S n • 


, - d^v d^v- 
hgsJi s 1; 


(50) 


With equation (50) and the definition of H, equation (47) can be reduced to 


Q-Cv) 


meC 




U \ 


(51) 


where has been replaced by m^. The variable on the right-hand side 

of (51) is related to v by equation (38). Equations (46) and (51) are the 
detailed balancing relations for this photoionization encounter. 


When these expressions are introduced into the gain expression (39) and 
the result combined is with the loss expression (41) and then integrated over 
all possible ion velocities. 



c^a 


j 


2Hv^p| 



Ss 



(52) 


With the aid of equations (40) and (46) we can rewrite (52) as 








ScoP 


sr £ 


d^v. 

1 


£ 


(53) 


This equation represents the net gain of electrons to the set per 

unit volume per unit time for the reaction represented by (37) . On summing 
over all possible species and states, we obtain the following collision 
integral for photoionization; 



n,£, j 



g 


S£*^£ 



(54) 
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IV. COLLISION INTEGRALS FOR NONELASTIC PROCESSES AFFECTING 
THE TRANSFER EQUATIONS OF THE HEAVY PARTICLES 
AND THE PHOTONS 


In an energetic plasma the heavy particles as well as the electrons are 
influenced by nonelastic collisions. Although this effect on the heavy parti- 
cles manifests itself primarily as a change in species concentration, momentum, 
and energy density rather than a change in the heavy particle distribution 
function itself, it will prove useful in the section on rate expressions to 
have at hand a convenient collection of nonelastic collision operators for the 
heavy particles as well as for the photons and the electrons. 

We will consider here radiative excitation and de-excitation reactions 
(indicated by subscript Ra) in addition to those reactions considered in the 
previous section. This reaction proves important in any calculation of number 
densities and energy transfers for the heavy particles and the photons and 
thus affects the electron distribution function in a coupled problem (ref. 9). 


The collision operators we will be considering in this section are 
indicated symbolically as follows: 

(1) Neutral species: 



(2) Ionic species:^ 



(3) Photons : 


- 



V " /ne - 

V /Ra 

\ /pj^ 


( 55 ) 


(56) 


(57) 


In this section, except for the development of the radiative excitation and 
de-excitation terms, we shall rely on the previous section for details and 
merely apply those results to the present formulation. 


^Collision integrals for ionic excitation or ionization can be deduced 
directly from the appropriate terms in equation (55) and will not be presented 
here . 
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Neutral Species 


As equation (55) indicates, the nonelastic collisions being considered 
that involve neutral species are: collisional excitation and de-excitation, 

collisional ionization and three-body recombination, radiative excitation and 
de-excitation, and photoionization and two-body recombination. All of these 
interactions except for radiative excitation and its inverse have been con- 
sidered in detail in the previous section. Hence by inspection or by follow- 
ing a development similar to that outlined in the previous section we can 
write the excitation and de-excitation collision operator for a neutral 
species in state j as 



Note that in this expression k is not restricted to being greater than j 
and the subscript s’ is used for the electron involved in the gain term in 
(58), Of course (58) could have been written with two collision operators on 
the right-hand side; one for k > j and one for k < j, for which s’ would 
be replaced by a 0 and a 1, respectively, consistent with equations (3) and 
(4) with subscripts nj and interchanged. 


We can also follow a development similar to that outlined in the previous 
section and write for the collisional ionization and three-body encounter: 


. ’3t 




3t 


Ion 


' Ion 


1 

H - f3f„.)g3 .o'. dfl23 d = 


(59) 


From Section III (p. 17) we can also deduce 



for the photoionization and two-body recombination encounters. 

Radiative excitation and de -excitation . - For the radiative excitation and 
de-excitation process we shall follow the development associated with photo- 
ionization and two-body recombination encounters (pp. 14-17), accounting for 
the differences between the ’’discrete” spectrum associated with this bound- 
bound transition and the continuous (above a limit) spectrum associated with 
the bound- free photoionization transition. Here, however, the development 
will be carried out with f^ rather than f^ for the photon distribution 
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function. A typical reaction for this case can be represented by 

"i k 

CAtom)^ + photon J fAtom)^ 

The energy equation for this encounter is 


( 61 ) 


2 hv = 2 mnVnj^ + 


(62) 


Let Qj (v) represent the cross section for the radiative excitation of a 
neutral atom from state j to state k as indicated in equations (61) and 
(62). The number of atoms in the set , d^Vj^. which are excited to 

atoms per unit time per unit volume as a result of radiative exci- 
tation can be represented by 


^nj d^Vnjf^(v,Xy)dv dr2^i|j^(v)cQ^ 


(63) 


In equation (63) (v) is a line-shape factor, which accounts for the fact 

that the frequency over which the absorption occurs is_^not precise but varies 
over a narrow range. The number of atoms in the set relax 

by radiative de-excitation to Vn*, d^v^. atoms per unit time per unit volume 
can be represented by 


d^Vnj^P^i(j^(v)dv dfi 


(64) 


As in expression (41) the differential cross section for relaxation 

(v, (v)dv d^2^ has spontaneous and stimulated components. Thus we can 

write 


Pj = aJ(v) + I(v,Xv)bJ(v) 


(65) 


where the first term on the right corresponds to spontaneous relaxation and 
the remaining term corresponds to the induced relaxation. The differential 
elements dv represent the range. and direction of the emitted photon. 

The line-shape factor for emission will be assumed to be the same as the 

^ ^ "1 k 

line-shape factor for absorption;^ thus w 

K 1 


Applying the principle of detailed balancing we can equate (63) with (64) 
and (65) at equilibrium and obtain 

Qf 




Ji 

'Eq hv 


( 66 ) 


^This assumption is strictly true for Local Thermodynamic Equilibrium 
(L.T.E.) and usually is assumed to hold for situations in which collisions con- 
tribute predominately to the excited state population distribution (ref. 18) . 
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In obtaining (66) equation (40) was used to replace by I in (63) . Also 

used was the differential volume element relation ^^^nj = which 

follows directly from a momentum balance and the fact that the momentum of the 
photons is neglected. Now using 


(I) 


^ 2hv^ 1 

Eq - 1) 


and the relations given by (10) and (11) of Section III (p. 3) we find 
from (66) 


A? = 

\ r2 k 


and 


k _ 

“j 2 v2 


4 


(67) 


( 68 ) 


The cross section Q. in (63) could be replaced by an Einstein coef fi- 
le ^ 

cient for absorption Bj such that (63) would be written as 

fji . d^VT^ . dv 

nj y 3 ^ 


This suggests the following relationships 

k Qi 


k 

A_ 

hv 


and 


the latter being well known. 


b" 


uUt 

(jOj k 


When these detailed balance expressions are introduced to the appropriate 
gain expressions and the result combined with the loss expression and inte- 
grated over the photon range and direction, one obtains 


or 


I 3t 


d^vn^ = d^V; 


/Ra 

j^k 


^9efnA 

— d^v_. = d^Vn. 
j-e->-k 



IB 








dfiv dv 


k-j 


(69) 


(70) 


21 



On carrying out the integration over the solid angle and summing over states 
we obtain 


or 


L~9t 




e"-nr 


/Ra 
k> j j 


3t 


/Ra 


j < k 



IB^ 


fn. ^ - fri JPv'f'v 

I "^3 pj ”k/ k^k 


Tt j 


/ e n • 


Ra 
j <k 


J 


3t 


k 

j <k 


/Ra 

j^k 




Note 


that by definition dv = 1 (ref. 17), 


(71) 


I Ph) dv (72) 


If we want to account for either excitation or de-excitation from a 
single state we combine (71) and (72) and write 





(73) 


For the second term on the right-hand side of (73) we would use (71) with 
k and j interchanged. 


Ionic Species 

The two interactions indicated by (56) have been considered in detail in 
Section III (pp. 8-17); thus by a slight intuitive extension of that work we 
can write directly: 


\~ir~ 













■S£ 


^3pS3Z 




IS 


d^v. 


(74) 


for the rate of change of the distribution function for ions of species n in 
state £ as a result of collisional ionization and three-body recombination, 
and 
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for the rate of change of f^ per unit volume per unit frequency as a result 
of radiative excitation and its inverse. 


It should be noted here that to the extent that the internal energy state 
of a particle is independent of its translational energy we can integrate over 
the heavy particle velocity in equation (78) directly and obtain 




(79) 


j 

k> j 


V. RATE EXPRESSIONS 


The rate of change of any electron property 4>(v5) as a result of non- 
elastic collisions can be designated as 




d^v = 


3e(nect-)' 


8t 


(80) 


Jne 


where (j> is the mean value of (j). Similarly, the rate of change of any 

property of the neutral species (|>(vn*), the ions or the photons 

3 ^ 

(|)^(v) as a result of nonelastic collisions can, respectively, be written as 




at 


'NE 






at 


(81) 


NE 


with corresponding expressions for the ions and the photons. When equa- 
tions (1), (55), (56), and (57) are combined with the above expressions. 


Ph 




3e(ne'f>)’ 


"aeC’^e'f’)' 

- 1 - 

-ae(nec^)' 

at 

NE 

_ at ^ 

Ex 


lon 

^ 9t 


(82) 


"SeCnn^nj)' 


^e (^nj ‘f’nj ) 

4 - 



ra (n 6 )1 

1 J 

4 - 

r9„(n_ cj) )~| 
3 3 

at 

NE 

at 

i 

Ex 

at 

Ion 

at 

Ra 

at 

— — 


"a (ni (|)^ )■ 



- 1 - 

^e C^ij^^ij^) 

at 

NE 

at 

lon 

at 


Ph 


(83) 

(84) 







at 

NE 

at 

Ra 

at 


Ph 


(85) 
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( 86 ) 


In regard to equations (83) and (84) we can sum over states and obtain 


and 


■3e(%<)>n)1 y 

^e ^nj ) 

L Jne 4- 

3t 


NE 


3 



(ni <t>i )1 
^n ^n 

-»NE 



(Hi^^- 

Jne 


(87) 


To check the consistency of the formulation of the collision integrals, in 
this section the terms on the right-hand side of equations (82)- (85) will be 
evaluated for various forms of the property functions. It will be demon- 
strated that equations (82)- (87) with the collision operators developed in the 
preceding sections will yield nonelastic expressions that agree with the usual 
conservation equations. In particular, for the reactions considered herein, 
if = 1 it will be shown that 


(- 


e^e 



0 


and 







( 88 ) 


The last two terms in the above equality are the rate of creation and destruc- 
tion of ions and atoms, respectively. Equations ( 88 ) state that excitation 
and de-excitation type collisions do not result in any net creation of elec- 
trons but that ionization and photoionization and their inverses are the only 
terms contributing to a net production of electrons. In addition, of course, 
the net rate of increase of free electrons equals the net rate of creation of 
ions which also equals the rate of destruction of neutral particles. It will 
also be demonstrated that if cj)(v 5 ) = (l/ 2 )m 0 V 3 ^, 



energy gained by neutrals, ions, and photons 
for the Ex, Ion, and Ph reactions 
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where (l/2)mgVg^ is the mean translator/ kinetic energy o£ the electrons. 

This expression shows that the net gain of thermal energy by the electrons 
during the encounters considered are a result of the following: 

(a) The net relaxation of excited states (that is, (9e’^nj/9t)gj^ j-f->-k 
the net production of state- j species from state-k species (k > j) and 

(b) The net production of neutral species, at the expense of the ions, 

(c) The net gain of photon energy, as a result of neutral species 
production, (3gnnj/9t)pj^^ and 

(d) The net loss of thermal and radiant energy, respectively, of the 
heavy species and photons involved in the encounters . 

The case when (fiCvg) = mgVg is not as interesting as the preceding cases 
and is briefly treated separately in the momentum transfer part of Section V. 


Inelastic and Superelastic Collisions 
Equation (19) can be used to write 

9g (ng^) 


at 


JEx 


-fds) d3v3 


gsj°sj ‘^^os d^Vg 


n, j,k 


(flS - fsfnk)«sk“s 


ffn, dfi, d^v„ d^v 

s nw®sk sk is s 


(90) 


For (}>(vs) =1 a typical term of equation (90) can be expanded: 


26 






fff ^ <»«„s ‘I’vj 

-ffJVnjSsi<^]) -30^5 d3v3 

*///*'"": ^ "“is "X ""''s 
-fIfVniSsV<^A ‘3“is ""''s 


(91) 


Combining the detailed balance relations, (13) and (14), and the relations 
between the differential velocity elements, (12) and (15), of Section III with 
this expression yields 

"“o= "'’'"j 
"“is "’’'l 

■j3J^s^ntgsk°k '*”ls ‘3®X '3^''s 

By this change of variables (v^^jVg Vj^^^,Vq and Vnj^,Vg ^ in the first 

and third terms of (91), respectively, the threshold energy dependence of the 
first set of integrals on the right-hand side of (91) has been removed and a 
threshold energy dependent cross section has been introduced into the third 
set of integrals of this equation. Now, since the subscripts associated with 
the electrons can be changed without affecting the value of the integrals, the 
first and second terms in equation (92) will cancel with the fourth and third 
tenns, respectively. That is, if in the first and second terms of equa- 
tion (92) we consistently change s -> 1, 0 ^ s and use the fact that 
d^si = will recreate the negative of the fourth and third 
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integrals.^ Thus, when <}>CVs) = 1, 


= 0 which implies : 


Ex, 



( 93 ) 


This is as expected since excitation and de-excitation encounters do not 
create or destroy electrons. 

For <t>(Vs) = (l/2)mgVs^ a typical term of (90) can be expanded by 
making use of detailed balancing expressions, energy equations, and variable 
transformations to obtain 


-(t 


at 


Ex 

j-^k 


“ ///& * 7 %(''hk 

- fffi "eVs"fsfnjgsj»sj 

- JJfl 


- + A., f„f g , o\ dfi 
^1/ iKJ ° nj^*ok ok os 

dJ^ d^Vn d^Vc 
os "j s 

- ''nj - ‘jkJfifn.Siyd 

dfi, d^v„ d^v_ 

Is u s 


d3y„ d^v 
"k o 


d^Vj, d^v, 

j 


(94) 


In this equation, in a manner analogous to that leading to equation (93), 
the first term in the bracket of the first set of integrals cancels with the 
last set of integrals, and the second set of integrals cancels with the first 
term in the bracket of the third set of integrals. Thus, 


change of variables of integration to a center of mass and relative 
velocity dependence will explicitly illustrate this equality. 
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{j HenieVe^^ 


dt 


=: A 


I Ex 
j^-^k 
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o^%^ok^ok 


d^Vn. 

os ^^k ^ 


- jOJfl^njgija^j df^is d^v,,^ 


^ffh ™n%fofnj,gok^Jk 

d^Vi 

■^fffj”'n^n.^lfn.gijo]j dS2^^ d^v^ 

-fffj 


Again if the variables and dummy subscripts are changed and the results of 
detailed balancing are used this expression can be written as 



^e^e^e^ 


at 


JEx 

j^k 




j^k 



( 96 ) 


where here we recall that k > j. The first term on the right-hand side of 
( 96 ) is obtained from equations ( 58 ), ( 81 ), and ( 83 ) when 4^n- ” ^ ^ two- 

state atom (say j and k) . It is easy to show, via detailed balance relation- 
ships and dummy variable changes, that for this particular j k type 
excitation and de-excitation reaction 



j^k 



j^k 


( 97 ) 
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as expected. The bracketed terms on the right-hand side of equation (96) are 
obtained from (58), when k > j, and the last four terms in equation (95). If 
now the states and species are summed and equations (58) and (86) are used, 



Equation (98) is interpreted as follows: The net gain of kinetic energy by 

the electron gas as a result of inelastic-superelastic encounters is equal to 
the gain of energy from a net de-excitation of excited species minus the net 
gain of kinetic energy by the neutral particles involved. 

The last term on the right-hand side of (98) is obtained via the 
following sequence of equalities: 
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where 


On - ^ 


n>N - 


"N n " ^ 


and 


1 V 


VxT — ^n^n^-n 

N ^ n n n 


Ionization and Three-Body Recombination Encounters 

The rate of change of an electron property ^(v^) as a result of 
collisional ionization and three-body recombination encounters can be 
expressed as 



d^Vs 


■ r [/#«•> ■ ^s^nj)gsj^sj 


n,«,, j 


(99) 

Before equation (99) is evaluated for various values of 4>(Vs), it will prove 
helpful to have at hand expressions for the collisional rate of change of num- 
ber density of the heavy particles involved as a result of this type of colli- 
sion. If we set (|)n. = “ ^n ~ ^in ” ^ equations (83) -(87) (and make 

use of (55), (56), (59), and (74)), we find the following expected results: 


(-P) -(-p) 


( 100 ) 


Ion 


Ion 
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That is, the net rate o£ creation of ions is equal to the net rate of destruc- 
tion of neutral atoms. These results follow directly once the variables are 
changed in equation (59) by the following detailed balance result 


Hgsj<^sj dVn = g2£g4£°2 


45, 






d3 


Vt 


and the dummy subscripts s, 2, 4 are changed in the resulting equation to 

1, s, and 3, respectively, and use is made of the fact that dfi, = dJ2„i. 

is ^ i 

For <P(Vs) =1 it is fairly easy to show by equations (59), (74), 

(81) (eq. (81) applied to the ions as well as to the neutrals), (55), and (56) 
that from a typical term of (99) we obtain 



j-^5. 








Equation (102) in combination with (100) yields the expected 







j^5, 


Similarly, on summing over all states and species, we find 



( 102 ) 


(103) 


(104) 


Equations (103) and (104) illustrate that for this type of encounter, the rate 
of creation of electrons equals the rate of creation of ions or the rate of 
destruction of neutral particles. 

For (|)(vs) = (l/2)meVs^ a typical term of equation (99) can be 
expanded by making use of equations (21), (22), (31), and (32) to obtain 



9e(j nemeVg^ Jj 


at 


Jlon 


me(v2^ + V4^) + Ajg, + 


2 ■ 2 


<^^2s 


■®[i 

^jnFi’ ■ ^3^^ •" 7 ™nVnj - ^-jSL ~ ^l^n^gu o[j dn^g <i\ d^j 


d'v3 d3v3 


(105) 


Close inspection of equation (105) when the dummy character of the electron 
subscripts are accounted for reveals that the sets of integrals containing 
(l/ 2 )me(v 2 ^ + V 4 ^) and -2 [ (l/2)meVs^] are equal in absolute value and 
therefore cancel each other. 

Using the following equality 
1 _ ... 2 .. 2-. . 1 - -2 


2 "'efv^ - Vg^) + 2 "'nVn^ " 


mivf 
1 ^SL 




= jrn^Cv^^ - Vg^) + jm^v^. - A., 


1 

-y- + 2 '"ev 


in the third set of integrals in (105) and recalling that the ionization 
cross section is symmetric in the subscripts s and 3 we find that the 
remaining terms in (105) containing electron kinetic energies cancel each 
other. We are finally left with 
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■ JJJJ d\ d^ViJ 

C106) 

If the electron subscripts in the first set of integrals in the first of 
the above bracketed terms are changed from 2, 4^ and s to s, 3, and 1, 
respectively, and the detailed balancing result given by equation (32) is used 
along with - dQ^^, this bracketed term becomes equal to 

- (3eni^/9 1) Jon j < ? ” C^e^e/^^) jon ^ similar manner the second 

1_ _1 ^ .1 - •__ -f__ _1_ T frv _ 2 ~l/i-\- | 


bracketed term in (106) can be shown to equal “| Ion y 

The third brack eted term in equation (106) can be shown to equal 
- jse [^(l/2)nnjmj^Vnj J/3t I when the subscripts in the first part of this 

bracketed term are changed from 1, s, 3tos, 2, 4 and the appropriate vari- 
able changes through equation (31) are made. Thus equation (106) can be 
written as 


^e^2“ ^e^e^e^ 
9F^ 


Pe^e 
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(107) 
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On summing over states and species and using (103) we obtain 


TigmeVe^ J 


3t 


Ion 



j,l,n 


3e(j 1 ^) 


8t 


Ion 




3t 


( 108 ) 


“^lon 


where 


and 


nj = 


E -i 


m. 



^i™i 


nimi 


E 

n 


niniiVi^ 


The sum over n here is interpreted as a sum over the ionic species of type 
n. Equation (108) illustrates that the net gain of kinetic energy by the 
electron gas as a result of collisional ionization and three-body recombina- 
tion encounters is equal to the net gain of ionizational energy due to recom- 
bination minus the net gain of kinetic energy by the heavy particles involved 


Photoionization and Two-Body Recombination Encounters 

As a result of photoionization and two-body recombination encounters, 
changes in an electron property <t)(v 5 ) can be accounted for on a rate basis by 
the following equation: 
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( 109 ) 
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Again it is convenient to have expressions at hand for the rate of change of 
number density of the heavy species and the photons involved in this reaction, 
From equation (60), (75) -(77), and (81) (eq. (81) applied to ions and photons 
as well as to neutrals), it is immediately apparent that when (|)^. = 1 

j ^ 

the following results are obtained: 


/^e^nj 

l'"TE 




^ Ph 





which obviously extends to: 



( 110 ) 


(111) 


That is, photons and neutrals are created at the same rate by radiative 
capture. 

For cf>(Vs) = 1 a typical term of equation (109) when compared with (75), 
(81) (eq. (81) applied to ions) and (110), will yield 


3t 


3t 


^Ph 





which extends to 





( 112 ) 
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Ph 


^e^i 

"Tt“ 


Ph 


^e^n 
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Ph 
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This expression can be summed over states to yield equalities similar to 
equation (104) . 


36 



in the form of the energy equation for this 


For 4 )(v 3) = (l/2)mgV3^ 
encounter, that is. 


1 2 1 2 1 2 . 
^ meVg^ = hv + y mnVj^. - y 


a typical term of (109) yields 


3e(“ neHieV^^) 
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Ph 
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-r - fsfiJgs^Pl ^ ‘I 
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d^v,- „ d^Vj. 


- A. 


^3 n ' 
e"e 
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When equation (50) is substituted into the first set of integrals in equa- 
tion (114) and the result is compared with (77) and (81) (eq. (81) applied to 
photons), we see that this set of integrals is equal to 


9g(nPhv) 

at 


-• Ph 





(115) 


which is defined as the negative of the rate of change of radiant energy per 
unit volume, u^, resulting from this interaction. _The radiant energy density 
of the photon gas is defined by (ref. 17) u^ = y*hvfP dv. 

Similar comparisons of equations (50), (59), (75), and (81) (eq. (81) 
applied to ions and neutrals), with the second and third sets of integrals of 
equation (114), reveal these integral sets to be 

r^e(^ rijijinnVnj ^ 

3t 


J Ph 
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and 


3e(i 


Thus equation Cl. 14) reduces to 


3e(i ^emeV) 


|V\ . 

^e(2 ) 

3t 

Ph 

( at i 

^ ' Ph 

dt 




‘j* (^) 

^Ph 

This expression can be interpreted for this interaction as: The net rate of 

loss of electron kinetic energy per unit volume is equal to the sum of net 
gain of energy by the photons, the neutrals, the ions plus the ionization 
potential times the net rate of production of electrons. This statement may 
appear to be misleading unless it is noted that when the net production rate 
of electrons is positive, the net rate of change of radiant energy will be 
negative . 

On summing over states and species we can write for equation (101) 
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Radiative Excitation and De-excitation 


The rate of change of a property of the neutral species as a result of 
radiative excitation and its inverse can be expressed as (see eqs. (86), (81), 
(71)-(73)): 




I 


From a typical term in equation ^ 119 ) or from equations (69) and (70), we can 
deduce 


since 


/^e^nj 

\^"^t 


^ Ra 
j-^k 



j^k 




This leads directly to the following simple conservation of species 
expression: 




j^k 


( 120 ) 


Applying these expressions to equation (119) yields 



j<k j“^k j>k 

On summing over all states j we now obtain 



( 121 ) 


( 122 ) 


That is, the total number density of a species does not change as a result of 
radiative excitation and de-excitation type encounters . In obtaining equa- 
tion (122) we made use of expressions like (120) 

for both j < k and k < j . (Since k and j range over the same number of 

states, there is always an even number of terms in (122) that will cancel in 

pairs according to equation (120).) 

If (f)^ = 1 in the first term on the right-hand side of equation (85) is 

combined with (78), and the result is compared with the first moment of equa- 

tions (71) and (72) , 
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which extends to 
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That is, the net rate at which photons are created is related to the net rate 
at which excited particles are destroyed. 

For a typical term in (118) can be written 
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For a two-state encounter, making use of the momentum balance mj^Vn- = i^n'^nv 
yields ^ ^ 


3e(nnjninVnj ) 
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On summing over all states j, making use of (124), the momentum balance for 
each specific encounter, and the relations used in obtaining (122) we conclude 
that the momentum of the neutral species does not change as a result of these 
nonelastic encounters; that is. 
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For d)j^^ = m^Vnj/2 we can use the energy equation (62) and write a 
typical term in equation (118) as 
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d^v. 




j>k 
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From equations (78), (81) (eq. (81) applied to photons), and (115) 


and 




j > k 




dv d^v- 




Also using equations (71), 


(72), (120), and d^Vn- = d^ 

J k 


we have 


and 


A.. 
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dv d^v 
11 n. 
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Substituting these expressions into equation (126) yields 
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Now on summing over states j the last two sets of sums on the right-hand 
side are symmetric and 



k,j j^k k,j k^j 

j <k j>k 


and 
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With these expressions, equation (131) when summed over the j states becomes 
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The first two sums on the right-hand side of (132) can, on comparison with 
(71) and (72), be written as 

£ 

k 

j<k 

and 
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These two terms can be combined, via equation (73), to obtain 
"1 
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This result can be combined with equation (132) to yield 
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(133) 


Equation (133) shows that the net gain of kinetic energy by the neutral 
species is equal to the net loss of energy by the photon gas and the net loss 
of internal energy from the upper excited states. 


Momentum Transfer 

If (|)(v 3 ) = m^Vg in the above rate expressions, and the appropriate 
momentum equation is used for each interaction, it is relatively easy to show 
that the momentum lost by the electrons in each type of nonelastic collision 
is gained by the heavy particles involved. We can summarize this result 
briefly as follows: 

Inelastic and superelastic encounters . - The momentum expressions for the 
reaction given by equation (2) which correspond to equations (3) and (4), 
respectively, are: 


meVs + m^v^j = m^v^ + m^v^^^ 


or 


meVi + m^v^^. = m^V3 + m^^v^^ 

When these equations are substituted into the rate equation (90) with 
^(Vs) = nigVs and appropriate variable changes are made. 
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where = (1/nij^) ^ This expression illustrates the momentum exchange 

n 

between the electrons and the neutral particles. 

Ionization and three-body recombination,- Here the relevant momentum 
equations are 


meVs + + m-v- ^ 

m^v-, + “ nioVc + m^v^ + mTV 
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When these expressions are combined with the momentum form of equation (99) 
and appropriate variable changes are made. 
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where Vj ^ (l/m^) ^ ^i^i- The electron has exchanged momentum with both the 


neutrals and the ions, 


Photo ionization and two-body recombination.- When the momentum of the 
photons is neglected, the momentum equation for this reaction is 


m^^Vj^. = m^Vi, + m^v 


3 


l^lc 




On substitution into the momentum form of equation (109) and changing 
variables there remains 


(^e^e^e) 
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Ph 


Again the electrons have exchanged momentum with both heavy species. 


Nonelastic Terms in the Equations of Motion 

Thus far in this section the rate expressions associated with various 
nonelastic collision operators have been developed. Here the preceding 
results will be summarized with the following relations between nonelastic 
terms that might appear in the equations of motion of a plasma composed of 
electrons, ions, neutral species and photons: 
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Species conservation: 



Mass conservation (p^ Hging, etc.)* 


1 

Q. 

0 


V /NE 

“t /ne V « ^eJ 


Momentum conservation (neglect photon momentum): 
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The last term in equation (137) represents 
This sum can also be represented by 


(3enf^/8t) + (agTij^/at) pj^ . 

(used in eq. (138)) since 


(8enN/3t)g^^P^ 0. 


47 


I 



( 142 ) 



These expressions agree with the usual macroscopic species conservation 
equations (ref. 19) , 

If we sum over species such that 


n = Un + ne H- nj 


P = Pn + Pe ^ Pi 


PV 5 P^Vf^ + PjVj + PgVg 


and 


PV2 5 


PjVi' 


Pe^e 


2 + 


we find 





Equations (143) also agree with the usual macroscopic (ref, 19) conservation 
equations v;hen the right-hand side is identified as the nonelastic rate of 
change of the average internal energy of the gas mixture. 
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Radiative Energy Loss 


A common phenomenon accompanying nonelastic collisions in plasmas is the 
loss of radiative energy to the boundaries of the plasma. The amount of radi- 
ative energy lost to the plasma depends on many factors such as the plasma 
density, its optical characteristics, its geometrical characteristics and the 
type of transition. This energy loss mechanism affects, in a coupled fashion, 
the electron distribution function and the species number densities as well as 
the plasma thermal energy. 

Physically, in a plasma in which the electrons are redistributing their 
energy among the other species the following might be expected: 

(1) Energetic electrons collisionally excite or ionize heavy particles. 

(2) These heavy particles in turn either undergo the inverse to the 
preceding encounters or relax radiatively, emitting photons. (The method of 
relaxation depends upon the relaxation rates for each process. This will be 
discussed in the next section.) 

(3) The resulting photons are then either absorbed by or lost to the 
plasma, depending upon the opacity of the system to the frequency of the 
emitted radiation. 

A loss of photons manifests itself in at least two significant ways: 

(1) The excited states of the heavy particles, from which the radiative 
relaxation occurs, tend to become depopulated. 

(2) The electron gas looses energy equivalent to the volumetric radiant 
energy lost from the plasma. 

Analytically the expressions developed in the preceding and subsequent 
sections remain valid. However, in application one must solve the equation of 
radiative transfer, coupled with the appropriate equations of motion, to find 
the local rate at which photons and the radiant energy leave the plasma. 

The rate at which photons are lost from the plasma can be related to the 
rate at which excited states of the heavy particles are depopulated as can be 
seen from the radiative transfer equation (ref. 17) and the results of 
pages 35 to 45. In addition the energy moment of the radiative transfer equa- 
tion yields a relation between the temporal and volumetric changes in the 
photon energy density and the collisional term given by equation (142). 

Detailed accounting of radiative losses is a formidable task. To 
circumvent some of the difficulties Holstein (ref. 20) introduced and others 
(ref. 9) used local, frequency-dependent, energy loss factors. Holstein, in 
particular, evaluated radiation-escape factors for resonance radiation for a 
few geometric shapes. 

These local radiation-escape parameters are defined for each transition 
as unity minus the ratio of photon absorption to photon emission. Thus, the 
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difference between the emission and absorption terms appearing in the 
preceding expressions would be replaced by the product of the particular 
radiation-escape parameter and the appropriate emission term. The value of 
these parameters (between 0 and 1) determines the local opacity of the plasma 
to a particular frequency of radiation. 


Reaction Rate Coefficients 

The usual reaction rate coefficients for a reaction which proceeds in a 
forward Kf or backward direction can be obtained directly from the 

preceding results. For example, consider the collisional ionization-three- 
body recombination reaction described by equation (20) : The local net rate of 

change in electron number density as a result of this reaction is determined 
by an equation of the form 


Oe^e\ 

9t I ~ ^lon^e^N “ %ec^e 
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From equations (74) and (104) we can deduce that 
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A comparison of (144) and (145) yields 
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If desired the detailed balancing expression given by equation (32) can be 
used to obtain Kf and/or in terms of ionization cross sections rather 

than recombination cross sections. These expressions clearly show the 
dependence of the reaction rate coefficients upon the distribution functions 
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and the cross sections. In addition one can deduce that the macroscopic 
detailed balance result, general but 

does hold at equilibrium. 

In a like manner rate coefficients for the other reactions considered 
herein can be obtained. Such coefficients for unidirectional reactions are 
particularly useful when competing reactions have dominant terms. In these 
cases one needs to consider only the rate coefficients for the dominant 
processes rather than the complete nonelastic collision integrals. 


VI. CONCLUDING REMARKS 


Collision integrals have been developed for some nonelastic processes 
that are important in plasma dynamic problems. On taking moments of these 
collision integrals, we have also demonstrated compatibility between these 
results and the nonelastic terms of the macroscopic equations of motion. 

The collision integrals developed here may prove to be unwieldy in their 
present form for many plasma dynamic analyses. Many simplifications are 
possible to make the collision integrals more tractable. The particular sim- 
plification employed depends upon the problem at hand. A few of the more 
successful simplifications, that is, the small electron mass approximation, 
the spherical harmonic expansion of the electron distribution function, the 
use of Maxwellian distribution functions, and the use of isotropic cross 
sections for nonelastic processes, are treated in some detail in the first 
16 references. 

It is hoped that the development presented herein allows a clearer 
understanding of the physics implicit in many less detailed treatments of 
nonelastic collision integrals. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, April 7, 1969 
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VII. APPENDIX 


NOMENCLATURE 


This is a partial alphabetical listing of the present nomenclature, 
symbols are defined locally in the text where they are first used. 


c 

e 

^m 

#(v) j 
fh 

®mh 

H 


coefficients for radiative de-excitation and excitation 
(eqs. (65) -(69)) 

speed of light 
electron 

electron velocity distribution function, f(^'^Tn) 
photon distribution functions (eq. (39)) 

velocity distribution function for heavy particle h 
relative speed [v^ - Vh| 

h^ “i 

parameter associated with Saha equation — r — ^ 


h 

T(v) 

|J| 

^Ion» ^Rec 
k 

mh 


’^e>’^h 


'■itSsP 


Planck’s constant 

specific intensity of v-radiation integrated over all solid 
angles 

Jacobian of a transformation, usually defined locally 

reaction rate coefficients defined on page 50 

Boltzmann’s constant 

electron mass 

mass of heavy particle 

number densities: electrons, heavy particles 

differential cross section for radiative capture (eqs. (41) 
and (42)) 


All 
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pj(v) 



T 

uR 


V ,V, 
h 


^iran* 




\ 


> 


/ 


3 


coefficient of differential cross section for radiative 
de-excitation (eqs . (64) and (65)) 

distribution functions associated with photionization and radiative 
excitation cross section (eqs. (39) and (63)) 

equilibrium temperature 
radiant energy density 

velocity of electrons, velocity of heavy particles 
electron velocity defined by (^m>Vm’^ 

coefficients for spontaneous and induced radiative capture (eq. (42)) 


parameter defined by, 


A 





excitation or ionization potential associated with a particular 
transition (eqs. (3), (21), (38), and (62)) 


6^ ^ delta function defined on pages 8-13 

memi 

U4 reduced mass, 

^ ^ m^ + mi 


V 


radiation frequency 


P 3 P e > 
Pn’^I 


mass densities defined on pages 46-48 (p^ = n^m^, etc.) 


k 3 
o 

mj mk 


angular distribution functions associated with inelastic and 
superelastic differential cross sections (eqs. (5) and (6)) 


h* 

^mh 

h* 

^mm*h 




distribution functions associated with ionization and three-body 
recombination cross sections (eqs. (23) and (24)) 


property of electron gas, <Kv,j) 
property of hea-vy particles, <|)(vj^) 
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R 


property of photon gas, <t>^(v) 

X^f angle of "deflection" of electrons 

X angle between relative velocity vectors g and g 

gs®3 sZ 3SL 

k i 

line shape factors for absorption and emission 

3 ^ 

Q solid angle used in describing particle direction 

degeneracies associated with states of an atom and its 

ion (o)^) 

Superscripts 

averaged quantity; velocity averaged or averaged over angles 
see subscripts j,k,£ 


C ) 

j,k,£ 


m 

m* 

n 

nj,nk 

i 

h 

I 

CM 

Eq 

NE 


Subscripts 

association with electron of speed Vj^^ (m = s, 0, 1, 2, 3, 4) 
association with electron of speed Vjj^’ 
neutral species of type n 

neutral particles of type n in states j, k 

ion associated with type n neutral 

ion in state £ associated with type n neutral 

heavy particle, i, n, n j , nj^, i^, I, N, and states j, k, £ 

ion mixture 

center-of-mass coordinates 
equilibrium 

quantities associated with nonelastic encounters 


Ex, Ion, 
Ra,Ph 


association with collisional excitation and de-excitation encounters, 
collisional ionization and three-body recombination encounters, 
radiative excitation and de-excitation encounters, and photo- 
ionization and two-body recombination encounters 


54 



1 


( ) 


( ) 


( ) 


Ex 

Ion 

Ra 


particular nonelastic encounter and its inverse which results in the 
heavy particle changing between a neutral in state j to a 
neutral in state k or an ion in state I 


C ) 


Ph 


j ,k,«. 


states of a neutral particle (j,k) and its ion 


Miscellaneous 


(Atom)^ 


atom of type n in state j 


^e 

3t 


rate of change of quantity per unit volume in phase space as a 
result of collisions 




volume elements in velocity space 
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TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons, 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 

Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non -aerospace 
applications. Publications include Tech Briefs, 
lYchnology Utilization Reports and Notes, 
and Technology Surveys. 



